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A  two-dimensional,  two-phase,  non-isothermal  model  was  developed  for  an  active,  tubular,  liquid-feed 
direct  methanol  fuel  cell  (DMFC).  The  liquid-gas,  two-phase  mass  transport  in  the  porous  anode  and 
cathode  was  formulated  based  on  the  multi-fluid  approach  in  the  porous  media.  The  two-phase  mass 
transport  in  the  anode  and  cathode  channels  was  modeled  using  the  drift-flux  and  the  homogeneous 
mist-flow  models,  respectively.  Water  and  methanol  crossovers  through  the  membrane  were  considered 
due  to  the  effects  of  diffusion,  electro-osmotic  drag,  and  convection.  The  model  enabled  a  numerical 
investigation  of  the  effects  of  various  operating  parameters,  such  as  current  density,  methanol  flow  rate, 
and  oxygen  flow  rate,  on  the  mass  and  heat  transport  characteristics  in  the  tubular  DMFC.  It  was  shown 
that  by  choosing  a  proper  tube  radius  and  distance  between  the  adjacent  cells,  a  tubular  DMFC  stack  can 
achieve  a  much  higher  energy  density  compared  to  its  planar  counterpart.  The  results  also  showed  that  a 
large  anode  flow  rate  is  needed  in  order  to  avoid  severe  blockage  of  liquid  methanol  to  the  anode  electrode 
due  to  the  gas  accumulation  in  the  channel.  Besides,  lowering  the  flow  rate  of  either  the  methanol  solution 
or  air  can  lead  to  a  temperature  increase  along  the  flow  channel.  The  methanol  and  water  crossovers  are 
nearly  independent  of  the  methanol  flow  rate  and  the  air  flow  rate. 
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1.  Introduction 

In  the  past  decade,  the  direct  methanol  fuel  cell  (DMFC)  has 
been  extensively  investigated  as  a  promising  power  source  can¬ 
didate  for  electronic  devices,  such  as  laptops  and  portable  digital 
assistants.  The  unique  advantages  of  the  DMFC  include  high  energy 
density,  low  emission,  instant  recharging,  and  compact  structure. 
However,  the  commercialization  of  the  DMFC  is  still  hindered  by 
several  technical  problems,  including  the  low  electro-activity  of 
methanol  oxidation  at  the  anode  as  well  as  methanol  and  water 
crossovers  through  the  membrane  [1-11].  Great  efforts  continue 
to  be  made  to  resolve  these  problems. 

Previous  DMFC  research  and  development  were  mainly  based 
on  the  planar  cell  structure  [2-14],  and  relatively  little  work  has 
been  reported  about  the  tubular  DMFC  [  1 5-24],  which  is  all  exper¬ 
imental  work.  Compared  to  the  planar  cell  structure,  the  tubular 
cell  structure  has  some  specific  advantages,  including: 
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(i)  it  has  a  larger  electrode  surface-to-volume  ratio,  and  thus,  a 
higher  system  energy  density  can  be  achieved; 

(ii)  it  has  a  greater  flexibility  in  shape; 

(iii)  it  eliminates  the  need  for  additional  plates  on  both  the  anode 
and  the  cathode  sides  to  hold  the  structure  together. 

Due  to  these  advantages,  the  tubular  cell  structure  has  been 
widely  used  in  solid  oxide  fuel  cells  (SOFCs).  However,  the  appli¬ 
cation  of  the  tubular  cell  structure  on  low  temperature  fuel  cells, 
such  as  polymer  electrolyte  fuel  cells  (PEFCs)  and  DMFCs,  has  rarely 
been  reported. 

Recently,  some  experimental  work  has  been  carried  out  to 
develop  tubular  DMFCs,  with  the  efforts  focusing  mainly  on  the 
methods  of  fabricating  the  tubular  membrane-electrode-assembly 
(MEA)  and  on  the  improvement  of  the  cell  performance.  Kunimatsu 
and  Okada  [16]  reported  the  performance  of  a  tubular  DMFC  using  a 
micro-tubular  electrolyte.  It  was  shown  that  both  the  catalyst  layer 
fabrication  and  hot-pressing  process  were  crucial  to  the  tubular  cell 
performance.  A  peak  power  density  of  12  mWcnrr2  was  achieved 
with  an  air-breathing  tubular  DMFC  fed  with  1  M  methanol  solu¬ 
tion,  which  indicated  that  in  terms  of  power  density,  tubular  DMFCs 
can  perform  as  well  as  planar  DMFCs.  Qiao  et  al.  [17-19]  reported 
different  methods  for  fabricating  the  tubular  anode  and  cathode, 
including  the  impregnation-reduction  method  with  heat  treatment 
for  the  anode  fabrication,  impregnation-reduction  method  with  a 
new  plating  process  for  the  cathode  fabrication,  and  a  wet  chemi¬ 
cal  assembling  process  for  the  formation  of  the  anode.  It  was  found 
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Nomenclature 

Alg 

interfacial  specific  area  between  liquid  and  gas 
phase,  m2  m-3 

A, 

specific  area,  m2  m-3 

Gw 

water  vapor  activity 

C 

molar  concentration,  mol  m-3 

Cp 

specific  heat  capacity,  J  kg-1  T-1 

D 

diffusivity,  m2  s-1 

F 

Faraday  constant,  96,478  C  mol-1 

Gr 

Grashof  number 

h 

heat  transfer  coefficient,  Wm-1I<-1;  enthalpy, 
Jmol-1 

frig 

interfacial  transfer  rate  constant  for  methanol, 
m2  s-1 

frm 

mass  transfer  coefficient,  m-2  s-1 

I 

current  density,  Am-2 

iP 

parasitic  current  resulting  from  methanol  crossover, 
Am-2 

J 

molar  flux,  mol  m-2  s-1 

Jo 

exchange  current  density,  A  m-2 

Ja 

anode  current  density,  Am-3 

jc 

cathode  current  density,  Am-3 

I< 

permeability  of  porous  material,  m2 

k 

thermal  conductivity,  W  m- 1  K-1 

fee 

condensation  rate,  mol  (atm  s  m3 )_1 

K 

evaporation  rate,  (atms)-1 

Kb  Kg 

relative  permeability 

m 

source  term  in  mass  conservation  equation, 
kgm-3  s-1 

m 

mass  flux,  kgm-2  s-1 

M 

molecular  weight,  kg  mol-1 

N 

mol  flux,  mol  m-2  s-1 

nd 

electro-osmotic  drag  coefficient 

Pc 

capillary  pressure,  Pa 

Pg 

gas  phase  pressure,  Pa 

Pi 

liquid  phase  pressure,  Pa 

Q 

heat  flux,  W  m-2 

R 

gas  constant,  J  (mol  K)-1 

R 

source  term  in  species  conservation  equation, 
mol  m  3  s-1 

R 

interfacial  species  transfer  rate,  mol  m-3  s-1 

^contact 

ohmic  contact  resistance,  £2  m2 

5 

source  term;  entropy 

s 

liquid  saturation;  entropy 

T 

temperature,  K 

U,  V 

velocity,  ms-1 

Vo 

thermodynamic  equilibrium  voltage,  V 

Vcell 

cell  voltage,  V 

X 

coordinate,  m 

r 

coordinate,  m 

Greek 

aa 

anode  transfer  coefficient  at  anode 

ac 

cathode  transfer  coefficient  at  cathode 

y 

reaction  order  of  ORR 

8 

thickness,  m 

£ 

porosity  of  porous  medium 

P 

overpotential,  V 

X 

water  content 

[l 

viscosity,  kgm-1  s-1 

p 

density,  kg  m-3 

o 

interfacial  tension,  Nnrr1;  proton  conductivity, 
m-1 

Superscripts 

eff  effective  value 

ref  reference  value 

sat  saturated  value 

*  in  equilibrium 

Subscripts 
a  anode 

c  cathode,  or  capillary 

cr  crossover 

dry  dry  membrane 

e  electrolyte,  or  evaporation 

g  gas  phase 

1  liquid  phase 

mem  membrane 

ML  methanol 

MV  methanol  vapor 

p  parasitic 

rl  relative  value  for  liquid  phase 

rg  relative  value  for  gas  phase 

T  temperature 

W  water 

we  dissolved  water 

WV  water  vapor 


that  improving  the  fabrication  process  of  the  electrodes  on  the 
tubular  membrane  was  effective  to  improve  the  tubular  DMFC  per¬ 
formance.  Shao  et  al.  [20,21  ]  developed  a  new  method  to  fabricate 
the  tubular  DMFC.  A  tubular  titanium  mesh  was  used  to  support  the 
catalysts  and  collect  the  current  from  the  cell.  The  catalyst  layers 
were  fabricated  using  a  dip  coating  method— a  recast  Nation  mem¬ 
brane  was  coated  on  the  cathode  electrode.  A  dilute  sulfuric  acid 
solution  was  used  as  the  electrolyte  to  reduce  the  contact  resistance 
between  the  anode  and  the  membrane.  It  was  also  shown  that  the 
tubular  DMFC  performed  similar  to  the  conventional  planar  DMFC. 
Yu  et  al.  [23]  used  a  porous  silica  pipe  as  the  support  to  fabricate 
the  tubular  electrolyte  membrane,  on  which  they  further  investi¬ 
gated  the  cell  performance  of  the  tubular  DMFC  and  tubular  PEFC. 
Lee  et  al.  [24]  investigated  the  performance  of  a  tubular  DMFC,  in 
which  the  tubular  MEA  was  fabricated  by  curling  the  conventional 
MEA  into  a  tube. 

Like  in  the  planar  DMFC,  a  number  of  physicochemical  phenom¬ 
ena  take  place  in  the  tubular  DMFC,  including  species  transport, 
momentum  transport,  and  multiple  electrochemical  reactions,  of 
which  a  better  understanding  is  essential  to  optimize  the  tubu¬ 
lar  cell  structure  and  improve  the  cell  performance.  Extensive 
efforts  have  been  made  to  develop  two-phase  mass  transport  mod¬ 
els  for  planar  DMFCs  in  order  to  quantify  the  complex  transport 
phenomena  and  elucidate  the  mechanisms  of  different  transport 
processes  [11,12,25-31].  Wang  and  Wang  [12]  modeled  the  DMFC 
using  the  mixture  multiphase  flow  model,  in  which  the  species  in 
the  liquid  and  gas  phases  were  assumed  to  be  at  the  thermody¬ 
namic  equilibrium  condition.  Yang  and  Zhao  [11,26]  presented  a 
two-phase  mass  transport  model  for  the  DMFC,  which  was  devel¬ 
oped  based  on  the  multiphase  flow  theory  in  porous  media.  Xu 
et  al.  [25]  developed  a  one-dimensional,  isothermal,  two-phase 
mass  transport  model  for  the  DMFC.  The  respective  effects  of  three 
water  transport  mechanisms  through  the  membrane,  as  well  as 
the  effect  of  the  MEA  design  on  the  water  crossover  and  cell 
performance,  were  investigated.  Rice  and  Faghri  [30]  proposed  a 
transient,  two-dimensional,  multi-phase  model  for  a  DMFC  with 
a  passive  fuel  delivery  system.  The  passive  delivery  system  uti¬ 
lized  a  porous  medium  to  passively  deliver  methanol  to  the  fuel 
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cell  while  controlling  the  concentration  of  methanol  at  the  anode 
side  to  limit  the  amount  of  methanol  crossover.  Rice  and  Faghri 
[31]  further  developed  a  thermal  model  to  study  the  thermal  and 
start-up  characteristics  of  a  passive  DMFC  system.  Their  model  also 
included  continuous  and  discontinuous  phase  limitations,  as  well 
as  a  probabilistic  spread  of  the  porous  properties.  Recently,  Xu 
and  Faghri  [3]  investigated  the  water  transport  characteristics  in 
a  passive  liquid-feed  DMFC  using  a  two-dimensional,  two-phase, 
non-isothermal  model.  The  effects  of  various  operating  param¬ 
eters,  such  as  current  density,  methanol  concentration,  and  air 
humidity,  as  well  as  the  effect  of  the  cathode  hydrophobic  air  filter 
layer,  on  the  water  transport  and  cell  performance,  were  investi¬ 
gated. 

Flowever,  with  regard  to  tubular  PEFCs  and  DMFCs,  only 
one  investigation  about  the  modeling  of  tubular  PEFCs  was 
reported.  Al-Baghdadi  [32]  developed  a  three-dimensional,  non- 
isothermal,  mass  transport  model  of  a  tubular  PEFC,  which  was 
used  to  understand  the  many  physicochemical  interactions,  com¬ 
plex  electrochemical  reactions,  and  transport  phenomena  that 
cannot  be  studied  experimentally.  The  model  was  also  used  to 
study  the  effects  of  various  material  parameters  on  the  cell  per¬ 
formance.  To  the  authors’  knowledge,  no  numerical  work  in 
relation  to  the  tubular  DMFC  has  been  reported  in  the  litera¬ 
ture. 

The  objective  of  this  work  was  to  develop  a  two-dimensional, 
non-isothermal,  two-phase  mass  transport  model  for  a  tubular 
liquid-feed  DMFC.  The  liquid-gas,  two-phase  mass  transport  in  the 
porous  anode  and  cathode  was  formulated  based  on  the  multi-fluid 
model  in  the  porous  media,  and  the  two-phase  mass  transport 
in  the  anode  and  cathode  channels  was  modeled  using  the  drift- 
flux  and  the  homogeneous  mist-flow  models,  respectively.  Since 
the  micro-porous  layers  (MPLs)  play  an  important  role  in  fuel  and 
water  management,  this  model  also  includes  MPLs  between  the 
diffusion  and  catalyst  layers  of  both  the  anode  and  the  cathode. 
Firstly,  an  evaluation  of  the  advantage  of  a  tubular  DMFC  system 
in  terms  of  energy  density  compared  to  a  planar  system  was  pre¬ 
sented.  Then,  the  model  developments  of  the  tubular  DMFC  and  the 
numerical  results  obtained  were  introduced.  Both  the  mass  trans¬ 
port  characteristics  and  the  heat  behavior  of  the  tubular  DMFC,  with 
different  current  densities  and  fuel  and  oxygen  flow  rates,  were 
investigated. 


2.  Energy  density  evaluation  of  the  tubular  DMFC 
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Fig.  1.  Sectional  views  of  (a)  a  planar  fuel  cell  stack  and  (b)  a  tubular  fuel  cell  stack. 


Accordingly,  the  ratio  of  the  reactive  surface  area  of  the  tubular 
fuel  cell  stack  to  the  planar  fuel  cell  stack  is: 

7T  x  (2 R  +  TMEa)  x  L  x  (W  x  H/(2R  +  nR  +  2T]v[ea)2) 
WxLx(H/(TMea  +  Tbp)) 

_  7T  x  (Tmea  +  Tbp)  x  (2 R  +  TMEa)  ^ 

(2  R  +  nR+  2TMEa)2 

For  a  standard  planar  fuel  cell,  the  thicknesses  of  the  bipolar  plate 
and  MEA  are  about  3  mm  and  0.6  mm,  respectively.  Fig.  2  shows  the 
variation  in  the  ratio  of  the  reactive  surface  area,  which  indicates 
the  ratio  of  the  stack  energy  density  when  the  cell  performance  of 
a  unit  area  in  the  planar  and  the  tubular  stack  is  assumed  to  be  the 
same. 

It  is  seen  that  the  ratio  depends  on  both  the  tube  radius  and 
distance  between  the  adjacent  tubular  cells.  For  instance,  when 


The  energy  density  of  a  DMFC  stack  can  be  increased  by  increas¬ 
ing  the  reactive  surface  area  per  unit  volume.  For  a  tubular  DMFC 
stack,  the  reactive  surface  area  depends  greatly  on  the  geometry 
size  of  the  unit  tubular  cell,  particularly  the  radius  of  the  tube  chan¬ 
nel.  To  elucidate  the  benefit  of  a  tubular  DMFC  stack  versus  a  planar 
stack  with  regards  to  the  energy  density,  the  reactive  surface  areas 
per  unit  volume  were  evaluated  for  each. 

Consider  a  standard  planar  DMFC  stack,  as  shown  in  Fig.  1  a,  hav¬ 
ing  a  cross-sectional  area  of  W (width)  x  H  (height)  and  a  length  of  L. 
The  total  cell  number  can  be  approximated  by  H/(TMEa  +  TBp X  where 
?mea  is  the  thickness  of  the  MEA,  and  TBp  is  the  thickness  of  the 
bipolar  plate.  Thus,  the  reactive  surface  area  corresponding  to  the 
unit  stack  volume  is  IV  x  L  x  ( H/(TMEA  +  TBp)).  While  for  a  unit  tubu¬ 
lar  DMFC  stack,  which  has  the  same  volume  size  as  the  standard 
planar  stack,  the  total  cell  number  depends  on  both  the  tube  size 
and  distance  between  the  tubes.  Assuming  the  shortest  distance 
between  the  electrodes  of  two  adjacent  tubular  cells  is  n  times  the 
radius  of  the  inner  tubular  channel  (R),  as  shown  in  Fig.  1  b,  the  total 
cell  number  can  be  expressed  by  Wx  H/(2R  +  nR  +  2TMEA)2.  Thus, 
the  reactive  surface  area  corresponding  to  the  unit  stack  volume  is 
tz  x  (2R+Tmea)  x  (W  x  H/(2R  +  nR+ 2 TMEA )2 ). 


Fig.  2.  Variation  of  the  ratio  of  the  reactive  surface  area  of  a  tubular  fuel  cell  stack 
to  a  planar  fuel  cell  stack  which  have  the  same  volume  (assume:  rMEA  =  0.6  mm; 
7bp  =  3  mm). 
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Fig.  3.  (a)  Schematic  of  a  tubular  DMFC  cross-section  and  (b)  the  computational 
domain. 


n  =  1,  the  ratio  first  decreases  sharply  and  then  slowly  continues 
to  decrease  with  an  increase  in  the  tube  radius.  When  the  radius 
becomes  a  critical  value  of  2.0  mm,  the  ratio  is  1.0,  which  means 
that  the  reactive  surface  areas  are  the  same  for  both  the  tubular 
and  planar  stacks.  When  the  radius  is  larger  than  the  critical  value 
(2.0  mm),  the  reactive  surface  area  of  the  tubular  stack,  as  well 
as  the  energy  density,  is  less  than  that  of  the  planar  one.  With  a 
decrease  in  the  distance  between  adjacent  tubular  cells,  the  critical 
value  is  increased.  When  n  =  0,  which  means  the  adjacent  tubu¬ 
lar  cells  come  into  contact,  the  critical  value  increases  to  about 
4.7  mm.  Therefore,  to  achieve  a  higher  energy  density  when  using 
a  tubular  DMFC  stack,  both  the  tube  radius  and  distance  between 
the  adjacent  cells  should  be  considered. 


(iv)  Both  the  gas  and  liquid  phases  are  continuous  in  the  porous 
layers,  and  the  flows  of  gas  and  liquid  are  all  controlled  by 
Darcy’s  law. 

(v)  Only  water  and  methanol  are  considered  the  condensable 
species. 

(vi)  The  membrane  is  impermeable  to  both  gas  and  liquid,  and 
only  dissolved  methanol  and  water  can  crossover  through  the 
membrane. 

(vii)  The  liquid  and  gas  phase  temperatures  are  the  same. 

3.1.  Transport  multiple  governing  equations 

3.1.1.  Mass  transport  in  the  anode  porous  region 

The  anode  porous  region  includes  the  anode  DL,  anode  MPL,  and 
anode  CL.  Due  to  the  coexistence  of  hydrophilic  and  hydrophobic 
pores,  both  liquid  and  gas  phases  exist  in  the  porous  region.  The 
general  governing  equations  of  the  mass  and  momentum  conser¬ 
vation  corresponding  to  each  phase,  as  well  as  the  conservation  of 
species,  are  given  by: 

Mass  :  V  ■  [p,ui]  =  =  %a  (liquid  phase) 

(2) 


v  '  [Pg“g]  =  9(/^Us)  +  7  8(rgfs)  =  mg,a  (gas  phase)  (3) 

where  u  =  ui  +  vk  represents  the  superficial  velocity  vector  based 
on  the  total  cross-sectional  area  of  the  fluids  and  porous  medium. 
The  rate  of  mass  generation  is  represented  by  m. 

k 

Momentum:  U\  =  -I<  —  Vpj  a  (liquid  phase)  (4) 

Mi 

k 

ug  =  -I<-^Vpg,a  (gas  phase)  (5) 

Mg 

where  I<  is  the  intrinsic  permeability  of  the  porous  medium  and  /cr 
denotes  the  relative  permeability  of  the  phase. 

Substituting  Eqs.  (4)  and  (5)  into  Eqs.  (2)  and  (3),  respectively, 
yields  the  governing  equations  for  the  anode  liquid  pressure  (pla) 
and  gas  pressure  (pg,a): 


d  ( P\Kkv\  3pU\  1  3  f  PlI<krl  dpU\  _  . 

dx  \  i±x  dx  )  r  dr  [  piX  dr  )  u 


(6) 


3.  Model  formulation 

This  section  presents  a  comprehensive  steady  state,  two- 
dimensional,  two-phase,  non-isothermal  model  of  a  tubular¬ 
shaped  DMFC,  as  illustrated  in  Fig.  3.  The  tubular  DMFC  includes 
an  anode  flow  channel  (ACF)  inside  the  tubular  cell,  anode  diffu¬ 
sion  layer  (DL),  anode  micro-porous  layer  (MPL),  anode  catalyst 
layer  (CL),  polymer  electrolyte  membrane  (PEM),  cathode  CL,  cath¬ 
ode  MPL,  cathode  DL,  and  an  annular  cathode  flow  channel  (CFC) 
outside  of  the  tubular  cell.  The  following  presents  the  model  to 
formulate  the  mass  and  heat  transport  in  different  regions  of 
the  computational  domain.  The  following  main  assumptions  were 
made  in  the  model: 

(i)  The  flow  is  symmetrical  about  the  axis.  Therefore,  all  3/3 0 
terms  are  neglected  in  the  governing  equations,  and  the 
governing  equations  for  transport  within  the  MEA  become 
two-dimensional. 

(ii)  Two  phase  flow  in  the  anode  and  cathode  channels  are  based 
on  one-dimensional  drift-flux  flow  and  homogeneous  flow 
along  the  direction  of  the  channel,  respectively. 

(iii)  The  porous  layers  are  homogeneous  and  isotropic. 


3  (  PgKkrg  3pg,a  \  ( r  ^g^<rg  m 

dx  y  fig  dx  J  r  dr  y  fi\  dr  J  ~  §,a 


Species:  ^ML,a)+1j(n^ 

dx  r  dr 


9  ( rjeff  9CML,; 

“  dx  l  ML  a  dx 


9CML,a  i  p 
ML, a  )  +KML,a 


9(UgC|yiv.ga )  1  f)(n'gCMv.ga )  _  d  I  nett  9QviV,ga 

dx  +  7  dr  -  dx  \  MV-sa  dx 


(8) 


i  ^  9  |  rneff  9Qviv,ga\  p 
+  rdr  ( rDMV,ga  J  +RMV,ga 

3(ugCwv,ga)  ,  1  9(n'gCwV,ga)  _  3  /  neff  9Cwv,ga 
dx  +  7  dr  “  dx  \  dx 


(9) 


,  1  9  Aneff  9CWV,ga 
+  r  9r  I  rLW,ga  9r 


+  ^wv„ 


ga 


(10) 


where  Cmu,  Qviv.ga.  anc*  ^wv.ga  are  the  liquid  methanol  concentra¬ 
tion,  methanol  vapor  concentration,  and  water  vapor  concentration 
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in  the  anode  side,  respectively.  Dff  represents  the  effective  diffu¬ 
sion  coefficient  of  species  i  and  R2  denotes  the  mole  generation  rate 
of  species  i. 

3 A. 2.  Mass  transport  in  the  cathode  porous  region 

Similar  to  the  anode  porous  region,  both  liquid  and  gas  phases 
exist  in  the  cathode  porous  region  which  includes  the  cathode  DL, 
cathode  MPL,  and  cathode  CL.  The  general  governing  equations 
of  the  mass  and  momentum  conservation  corresponding  to  each 
phase,  as  well  as  the  conservation  of  species,  are  given  by: 

Mass  :  V  ■  [p,ui]  =  +  y  =  m,  c  (liquid  phase) 

(ID 


v  '  [Pg“g]  =  9(/^Us)  +  7  9(rg^s)  =  rhg, c  (gas  phase)  (12) 
k 

Momentum:  U\  =  -I< — Vpt  c  (liquid  phase)  (13) 

Mi 

Ug  =  -/C^Vpg,c  (gas  phase)  (14) 

Mg 

Similarly,  the  governing  equations  for  the  cathode  liquid  pres¬ 
sure  (plc)  and  gas  pressure  (pg,c)  can  be  obtained  from  Eqs. 
( 1 1 )-( 1 4),  which  give : 


d_ 

dx 

d_ 

dx 


( p\Kkr\  dp\c\  _  \_d_  f  PlI<krl  dplc\ 
y  Mi  9x  J  r  dr  y  fi\  dr  J 

f  Pg/^rg  dPg,c\  ( r  PgKkyg  9pg,c 

l  fig  dx  J  r  dr  l  fi\  dr 


=  mlc  (15) 

J  =  riig?c  (16) 


Species  : 


9(ugCo2,gc)  1  d(ri/gCo2,gc) 

dx  r  ^ 


Deff  ^O^gcA 

ax  ^  °2’gc  dx  J 


+ 


i  d_ 
r  dr 


dQ)2,gc 

dr 


+  ^02,gc 


(17) 


9(ugCwv,gc)  ,  1  9(rygCwv,gc)  _  9  /  neff  ^Cwv,gc  \ 

a7  r  a?  “  a7  ax  y 


+ 


1  _a 

r  9r 


eff 

WV,ga 


acwv,gc 

ar 


+  ^WV,gc 


(18) 


where  Co2)gc  and  Cwv.gc  are  the  gas  oxygen  concentration  and  water 
vapor  concentration  in  the  cathode  side,  respectively. 

Here  it  should  be  noted  that  the  capillary  pressure  for  the  porous 
region  is  given  by  [33]: 

Pc  =Pg  —  Pi  =o  cos  0(e/K)a5J(s)  (19) 

where  J(s)  represents  the  widely  used  Leverett  function: 

f  1.417(1  -s)-  2.120(1  -s)2  +  1.263(1  -  s)3  0  <  0  <  90° 

J(s)={  (20) 

1.417s  -  2.120s2  +  1.263s3  90°  <  0  <  180° 

When  the  gas  and  liquid  pressures  are  solved  in  every  iteration, 
the  liquid  saturation  is  determined  by  Eqs.  (19)  and  (20).  Recently, 
some  efforts  have  been  devoted  toward  the  determination  of  the 
capillary  property  of  the  fuel  cell  porous  media.  However,  there 
are  inconsistencies  in  the  newly  developed  capillary  functions,  and 
thus,  the  Leverett  function  is  still  used  in  this  work  [34]. 


3 A. 3.  Mass  transport  in  the  membrane 

Due  to  its  extremely  low  permeability,  the  membrane  is  usually 
regarded  as  a  gas  insulator.  However,  liquid  water  and  methanol 
in  the  porous  region  can  dissolve  into  the  electrolyte  with  the 
effect  of  adsorption  to  form  the  dissolved  phase  and  cause  trans¬ 
port  between  the  anode  and  cathode  through  the  membrane.  The 
transfer  of  the  dissolved  water  through  the  membrane  depends  on 
molecular  diffusion,  electro-osmotic  drag  and  convection.  Accord¬ 
ingly,  the  molar  flux  of  water  crossover  (NH2ocr)  through  the 
membrane  can  be  given  by: 

Nh2 Ocr  =  -Dwe(A)VCwe  +  nd, ^07  -  (21 ) 

Thus,  the  governing  equation  for  the  dissolved  water  concen¬ 
tration  (Cwe)  is: 

^  •  l^H2Ocr  —  V  •  (—  Dwe(A.)VCwe)  +  V  •  ^Ud, H20^  j  —  0  (22) 

Note  that  the  dissolved  water  concentration  Cwe  can  be  trans¬ 
formed  to  the  water  content  X  in  the  electrolyte  and  the  relationship 
between  Cwe  and  X  is  given  by  X  =  Cwe{EWl pdry). 

The  molar  flux  of  methanol  crossover  (NMLcr)  through  the  mem¬ 
brane,  which  also  depends  on  molecular  diffusion,  electro-osmotic 
drag  and  convection,  can  be  given  by: 


NMLcr  =  -^M,nVCML  +  ^d,M  ~F  ~  ^  ^1,Ca 


Ml<5n 


Cml 


(23) 


where  A p\  c_a  represents  the  liquid  pressure  difference  between 
the  cathode  and  the  anode. 


3 A  A.  Energy 

The  energy  equation  used  in  the  entire  computational  domain 
can  be  expressed  as  follows: 

d((Pgcp,gug  +  P\Cp,\U\)T)  1  d(r(PgCp,gVg  +  Picp,ivi)T) 

dx  r  dr 


d_ 

dx 


i  a  (  ieffar 
7S  <57 


+  Sj 


(24) 


where  kjf  represents  the  effective  thermal  conductivity  of  the  heat 
transfer  media,  and  Sj  denotes  the  heat  generation  rate. 


3 A. 5.  Anode  flow  channel  (inside  the  tube) 

Since  complicated  two-phase  flow  patterns  in  the  DMFC  anode 
and  cathode  channels  have  not  yet  been  established  quantitatively, 
a  one-dimensional  drift-flux  flow  model  and  a  homogeneous  mist 
flow  model,  along  the  flow  direction,  are  used  to  describe  the  two- 
phase  flow  in  the  anode  and  cathode  channels,  respectively.  The 
general  governing  equations  of  the  mass,  momentum,  species,  and 
energy  conservation  corresponding  to  each  phase  are  given  by: 


Mass  conservation  : 


9(A(1  ~  sg)ui) 

dx 


2^1, afc 

7j 


(25) 


d{PgSgUg)  _  2mg  afC 

dx  ~  r\ 


(26) 


where  sg  represents  the  cross-sectional  averaged  gas  void  fraction 
in  the  anode  channel,  U\  and  ug  denote  the  cross-sectional  aver¬ 
aged  velocities  of  the  liquid  and  gas  phases,  respectively,  and  mi  afC 
and  mg  afc  represent  the  liquid  and  gas  mass  fluxes  at  the  interface 
(r = ri )  toward  the  anode  channel. 


.  a  A  1  -  Sg  U,Ul  df  OaS  qUqUa ) 

Momentum  equation  :  — _  g  1  1  +  gg  -  — 

dx  dx 

=  -^-2pm-p?-um\um\ 

dx  Uh,nfc  '  1 


(27) 
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Table  1 

Boundary  conditions  for  the  tubular  DMFC  model. 


Location 

Variables 

Mass/molar  flux 

r=r  i 

Pll  +  =Pm,  Si|+  =(1  -Sg), 

mi\+  =mi,afC 

CmL  |  +  =  CMeOH , . 

mg  +  =  mg)afC 

T  |  +  =  <7afc/^Ta  +  Fafc 

l^MeOH  |  +  =  iVMeoH,afc 

N,|+=lVUfc 

r=r4 

II 

O 

*1? 

II 

O 

mil-  =0 

^we  |  +  — ^we,i>  ^We  1  ^we,v)^l- 

mg  |  _  =  o 

Ni\-  =  0 

r=r5 

CMeOH  |  +  =  0,  |  ^  =  0, 

m\+  =  o 

mg  =0 

■^we  |  =  ^we,v  (^we  1  ^‘we,v)^l  + 

N,|+=0 

r=r8 

Pg  —  Pm?  Sj |_  =  (1  —  Sg), 

mfl_  =  mi)CfC 

Co2,gc|  =  Co2,cfc?> 

mg  =  mg)Cfc 

T\-  =  Qcfc/tljc  +  fcfc 

Ni|_=NifCfc 

0<r<ri  and  (x  =  0) 

CMeOH  |  +  =  CMeOH, in  >  Pm  |  +  =  1  atm 

T  =  Tin,  fll  |  +  =  Ma,in>  %  |  +  —  0 

r8  <r<r9  and  (x  =  0) 

Co2  =  Quin,  Cwv  =  Cwv.in, 

Prn  |  +  =  1  atm,  T  =  Tin 

Ug|  =Uc,in,  Sl\+  =0 

r\  <r<r8  and  (x  =  0 
orx  =  L) 

+ 

II 

o 

+ 

II 

o 

mi|+  =0 

mg  |  +  —  o 

Ni|+=0 

^MeOH  |  +  =  0 

Heat  flux 


Q 1+  —  7afc 


Q  L  =  <7l+ 


Q  L  =  <71+ 


<7l_  =  Qcfc 


<71+  =  0 


where  pm  represents  the  pressure  inside  the  anode  channel,  and  um 
represents  the  volumetric  velocity  of  the  liquid-gas  mixture,  and 
is  defined  as: 

Um  =  0,(1  —  Sg)  +  OgSg  (28) 

The  relation  between  the  gas  velocity  and  liquid  velocity  is 
defined  by  the  drift- flux  model  [11,12]: 


3. 1.6.  Cathode  channel  (outside  the  tube) 

Similar  to  the  anode,  the  governing  equations  for  the  homoge¬ 
neous  mist  flow  in  the  cathode  channel  are  listed  as  follows: 


9(p,s,ui) 

Mass  conservation  :  * 

dx 


2  r8m1>cfc 

crI  -  rI) 


(33) 


9(P g(l  -S|)ug)  =  2 r8mg  ,cfc 

^  “  (r| -r|) 


(34) 


0g  =  C0Om  +  Ugj 


(29)  ,,  „  d{p\S\U\U\)  3(pg(l  -  s,)ugUg) 

v  ;  Momentum  equation  :  '  +  ^ 

3x  3x 


where  the  distribution  parameter,  C0  =  (1.2  -  0.2^/pg/p,)(l  - 
0-18. 0sg )^  an(j  t^e  drift  velocity,  u^  =  0.37 ^/2r^(p\  -  pg)g/p\,  are 
suitable  for  a  bubbly  flow  in  a  pipe  [35]. 


Species  : 

9(Sg0gC,,afc) 

9x 


9((1  -5g)0,CMeOH)  _  2NMe0H,afc 


dx 


2  Ni 


i,  afc 


(30) 


(31) 


where  Ni  afc  represents  the  molar  flux  of  the  species  (i.e.,  CML,a 
CMv,ga>  CWA)ga)  at  the  interface  (r= r\ )  toward  the  anode  channel. 


Energy  : 


d((Pgcp,gug  +  ACpjU,)Tafc)  _  2qx,afc 

dx  ~  ri 


(32) 


where  fafc  represents  the  cross-sectional  averaged  temperature  of 
the  mixture,  and  qTafC  represents  the  heat  flux  at  the  interface 
(r  =  ri )  toward  the  anode  channel. 


-  9^m  2  o 

~  2pn 


-Um  \  Un 


-^h,cfc 


(35) 


Species  : 


9((1  Si)ugQ?cfc)  _  2rsNi  Cfc 

^  “(r|-r|) 


Energy  : 


^((PgcP,gug  +  P\cp3u\)tcfc) 

dx 


2l~8QT,cfc 

(rl-rI) 


(36) 

(37) 


3.2.  Sub-models 

In  the  anode  of  the  DMFC,  the  kinetics  of  the  methanol  oxi¬ 
dation  reaction  (MOR)  is  modeled  by  the  widely  used  Tafel-like 
expression: 

ja  =  A/,aJoeMeOH  ( ^  eXP  (~pj (38) 
V  LMeOH  /  V  7 

where  the  reaction  order  (y)  is  related  to  the  methanol  con¬ 
centration  and  is  assumed  to  be  zero  order  when  the  methanol 
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Table  2 

Constitutive  relations  in  the  governing  equations  for  the  tubular  DMFC  model. 

Parameters  Expressions 


Relative 

permeabilities 

Effective  diffusion 
coefficients  of 
species 


Effective  thermal  conductivity  in  the  porous  regions 


Generation  rate  of 
mass  in  liquid 
phase 


Generation  rate  of 
mass  in  gas  phase 


Mole  generation 
rate  of  species 


Heat  generation  rate 


Table  3 

Cell  geometric  dimensions  and  operating  parameters  used  in  the  tubular  DMFC  model. 


Parameters 

Symbols 

Value 

Unit 

Tube  inner  radius 

r\ 

i.o  x  nr3 

m 

Anode  DL  thickness 

r2  —  IT 

2.6  xlO-4 

m 

Anode  MPL  thickness 

r3-r2 

0.3  x  10-4 

m 

Anode  CL  thickness 

r4-r3 

0.2  x  10-4 

m 

Membrane  thickness  ( Nation  115) 

r5-r4 

1.25  xl0“4 

m 

Cathode  CL  thickness 

r6-r5 

0.2  x  10-4 

m 

Cathode  MPL  thickness 

r7-r6 

0.3  x  10-4 

m 

Cathode  DL  thickness 

r8-r7 

2.6  xl0“4 

m 

Flow  channel  length 

L 

5.0  xlO-2 

m 

Annular  cathode  channel  size 

r9-r8 

1.0  x  10“3 

m 

Fuel/Air  inlet  temperature 

Fin 

333.15 

K 

Anode  channel  inlet  pressure 

Pi, in 

1.013  x  10“5 

Pa 

Cathode  channel  inlet  pressure 

Pg,in 

1.013  x  10-5 

Pa 

Inlet  methanol  concentration  at  anode 

Qvi.in 

1000 

mol  m3 

Inlet  oxygen  concentration  at  cathode 

Co2,in 

0.21  x  Pg/RT 

mol  m3 

Inlet  water  vapor  concentration  at  cathode 

Cwv.in 

0 

mol  ITT3 

Inlet  carbon  dioxide  concentration  at  cathode 

Cc02,in 

0 

mol  m3 

Inlet  gas  saturation  at  anode 

h 

0 

- 

Inlet  liquid  saturation  at  cathode 

Si 

0 

- 

kri=s4-5  Liquid 
/<rg  =  (l  -s)4-5  Gas 

Dfg  =  D,)ge1-5(1  -  s)1'5  i :  02,  WV,  MV 


Deff  = 


A 


5(£+£e) 


kf  = 


mu 


m,c 


[£/(DMl£l-5sl-5)+ee/(DMiN^-5)]DM!N 
esk\  +  e(l  -  s)kg  +  /<dl/mpl/cl 

^mem 

Mh20^v1  -  ^M^MV.ga 

M„2o  («v1  -  ^  -Mm  (rMv, 


(ADL(AMPL)/ACL/PEM) 

(DL(MPL,CL)/PEM) 


I  k  |  fP  ) 
ga  +  6F  +  6  FSX,  ) 


(ADL(AMPL)/ACL) 


Mh20^v1 


Mh2o  (  RV1  + 


[20  ^ 


Nh20ct  (jc-h/Scd) 
Sm  2  F 


W8cc\ ) 


(CDL(CMPL)/CCL) 


^02,gc  — 


^ML,a  = 


_  f  -Ru 

,a  —  S 

[  -Ru 


-Mh20^v1  +  Mivi^MV.ga 
-Mh20^v1  +  lWM^MV,ga  +  Mco2^C02,ga 
-Mh20^v1 

-Mo2Jc/ 4 F  +  Mco2  IP/6F8CC\  -  Mh2oRv\ 

,  Rwv,gc  =  ~Rv  l  (CDL(CMPL)/CCL) 


3  F8CC 
(ADL(AMPL)/ACL) 

(CDL(CMPL)/CCL) 


(ADL(AMPL)/ACL) 


-Jc/4F 
-FMV,ga 

-^MV,ga  -ja/6F 

Fwv,ga  =  -Fvi,  ^Mv,ga  =  FMv,ga  (ADL(AMPL)/ACL) 

Rv\Ahv  -  RMv,gAhMv  AD  L(  AM  PL) 

Ja(Pa  —  FZ\Smor/6F)  +  J2/ CTmem  +  Rvl^hp  —  R]VIV,g2l/lMv  ACL 
I^/cfmem  PEM 

jcOpd  -  FASorr/4F )  -  JpTASmor/6F<5cci  CCL 

+f2/  ^mem  +  RvlAhv 

<  RviAhv  CDL(CMPL) 


Sj  =  < 
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Table  4 

Physicochemical  properties  used  in  the  tubular  DMFC  model. 
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Parameters 

Symbols 

Value 

Unit 

Ref. 

Porosity,  permeability 

ADL 

£adl-  f^adl 

0.75, 1.0  x  10“12 

-,  m2 

[25] 

AMPL 

£ampl»  f^ampl 

0.3,  2.5  x  10-13 

-,  m2 

[25] 

ACL 

£-acl »  f^acl 

0.3, 1.5  xlO-14 

-,  m2 

[25] 

MEM 

£mem>  Kmem 

0.3,  2.0  x  10-18 

-,  m2 

[25] 

CCL 

£ccl>  Keel 

0.3, 1.5  xlO”14 

-,  m2 

[25] 

CMPL 

£cmpl>  Kcmpi 

0.3,  2.5  x  10-13 

-,  m2 

[25] 

CDL 

£cdl.  Kcdi 

0.75, 1.0  x  10“12 

-,  m2 

[25] 

Nation  volume  fraction 

ACL 

£e,acl 

0.3 

- 

[25] 

CCL 

£e,ccl 

0.3 

- 

[25] 

Diffusivities 

MeOH  in  water 

^M.l 

1.58  x  10-9eao2623(T-298) 

m2  s_1 

[25] 

MeOH  in  Nation 

Dm,N 

4.9  X  1 0 — 1 0 0l2436(1/333 — l/T)] 

m2  s_1 

[25] 

Methanol  vapor 

— 6.954  X  10-6  +  4.5986  x  10“8r+ 
9.4979  x  10_11T2 

m2  s-1 

[25] 

-  /  _  \  1.823 

02  in  gas 

f^02,gc 

1.775  x  10-5  (  27535  ) 

m2  s_1 

[25] 

Water  vapor 

Dwv.g 

2.56x10  ^  307T5  ) 

m2  s-1 

[25] 

Dissolved  water 
in  Nation 

13  we 

4.17  x  lO-^ieie-1  +  l)e-2436/T 

m2  s_1 

[25] 

Thermal  conductivity  of  membrane 

k  mem 

0.2 

Wm-1  K”1 

[29] 

Thermal  conductivity  of  DL  and  CL 

koi (kc  l) 

1.5 

Wm-1  K"1 

[29] 

Thermal  conductivity  of  MPL 

/<MPL 

1.5 

Wm-1  K"1 

[29] 

Thermal  conductivity  of  gas 

kg 

0.026 

Wm-1  K1 

- 

Thermal  conductivity  of  liquid 

kx 

0.62 

Wm-1  K-1 

- 

Heat  capacity  of  liquid  water 

Cpl 

4200 

Jkg-MC1 

- 

Heat  capacity  of  gas 

cpg 

1007 

Jkg-1  K_1 

- 

Viscosity  of  gas  phase 

Mg 

2.03  x  10“5 

kgm-1  s_1 

[11] 

Viscosity  of  liquid  phase 

Mi 

4.05  x  10“4 

kgm-1  s_1 

[11] 

Absolute  entropy  of  liquid  methanol  (1  atm,  298  K) 

s° 

MeOH,l 

126.8 

J  mol-1  K"1 

- 

Absolute  entropy  of  liquid  water  (1  atm,  298 1<) 

s° 

H20,1 

69.95 

J  mol-1  K"1 

- 

Absolute  entropy  of  C02  (1  atm,  298  K) 

s° 

C02 

213.685 

J  mol-1  K-1 

- 

Absolute  entropy  of  02  (1  atm,  298 1<) 

s9 

02 

205.033 

J  mol-1  K"1 

- 

Electro-osmotic  drag  coefficients  of  water  and  methanol 

nd,H20 

2. 5  \ 

22  A 

- 

[25] 

^d,M 

nd,H20*M 

- 

[25] 

Evaporation  rate  constant  for  water 

fee 

5 

(atms)-1 

Assumed 

Condensation  rate  constant  for  water 

kc 

5.0  x  104 

s-1 

Assumed 

Interfacial  transfer  rate  constant  for  methanol 

h\g 

0.05 

m2  s-1 

[11] 

Specific  interfacial  area  between  liquid  and  gas 

Aig 

105 

m-1 

[11] 

Proton  conductivity  in  membrane 

^mem 

7_3e[1268(l/298-l/T)] 

^2_1  m_1 

[25] 

Henry  law  constant  for  methanol 

k+M 

0.096e°°4511(T-273) 

atm 

[25] 

-2.1794  +  0.02953(T  -  273) 

The  saturation  pressure  of  water  vapor 

logioP^ 

-9.1837  x  10“5(T-273)2 
+1.4454  x  10~7(T  -  273)3 

atm 

[11] 

The  saturation  pressure  of  methanol  vapor 

nsat 

^MV 

k\\Xu\ 

atm 

[11] 

Latent  heat  of  methanol  evaporization 

A/lMv 

37.7  x  103 

J  mol-1 

- 

Latent  heat  of  water  evaporization 

Ahv 

44.9  xlO3 

J  mol-1 

- 

Thermodynamic  voltage 

Vo 

1.21 

V 

[25] 

Transfer  coefficient  of  anode 

Ua 

0.52 

- 

[25] 

Transfer  coefficient  of  cathode 

ac 

1.0 

- 

[25] 

Anode  exchange  current  density 

A  iref 

lxl03  exp  (^p  (£-*)) 

Am-3 

[25,29] 

Cathode  exchange  current  density 

A  iref 
/lv-cJ0,02 

2111  exp  (Z3200  (_i_  _  i)) 

Am-3 

[29] 

Anode  reference  concentration 

+ef 

100 

mol  m3 

[12] 

Cathode  reference  concentration 

36.5 

mol  m3 

[11] 

Surface  tension 

a 

0.0644 

N  m1 

[25] 

Equivalent  weight  of  ionomer 

EW 

1.1 

kg  mol-1 

[11] 

Dry  membrane  density 

Pdry 

1980 

kg  m-3 

[11] 

Contact  resistance 

^contact 

0.45 

£2  cm2 

Assumed 
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Fig.  4.  Distribution  of  methanol  concentration  in  the  anode  electrode  and  mem¬ 
brane  and  oxygen  concentration  in  the  cathode  electrode  (methanol  solution:  1  M; 
current  density:  120mA cm-2;  anode  flow  rate:  0.015 ms-1 ;  cathode  flow  rate: 
1.0ms-1). 


Fig.  6.  Distribution  of  liquid  saturation  in  the  porous  region  at  different  channel 
positions  (methanol  solution:  1  M;  current  density:  120  mA cm-2;  anode  flow  rate: 
0.01 5  ms-1;  cathode  flow  rate:  1 .0  m  s-1 ). 


concentration  is  higher  than  a  reference  value  [12].  Otherwise,  the 
first-order  reaction  is  applied. 

The  cell  current  density  can  be  calculated  by 

/  =  [  h  dx  (39) 

Jacl 

The  rate  of  methanol  crossover  is  expressed  by  the  ‘parasitic’ 
current  density: 


Jp  =  6FNMLcr 


(40) 
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Fig.  5.  Distribution  of  liquid  saturation  in  the  porous  region  (methanol  solution: 
1  M;  current  density:  120mA  cm-2;  anode  flowrate:  0.015  ms-1;  cathode  flowrate: 
1.0ms-1). 
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Fig.  7.  Distribution  of  cell  temperature  rise  (methanol  solution:  1  M;  current  den¬ 
sity:  120mA cm-2;  anode  flow  rate:  0.015ms-1;  cathode  flow  rate:  1.0ms-1). 
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where  is  the  saturation  pressure  of  water  vapor,  and  ywv  is  the 
molar  fraction  of  water  vapor  in  the  gas  phase.  The  mass-transfer 
coefficient  (ftvi) can  be  given  by: 

h  _  kce(l -s)ywy  A  |  |ywvPg-P^|\ 

'  !  2RT  ywvPg-p^,  J 

fceespi  /  |ywvPg-p^|\ 

+  2Mh2o  ^  ywvPg-P^,  j  1 

where  kc  and  ke  are  the  condensation  and  evaporation  rate  con¬ 
stants. 

The  rate  of  methanol  condensation  and  evaporation  between 
the  liquid  phase  and  vapor  phase  is  modeled  by  the  following 
expression: 

%v,g  =  AgM1  - s)(P%tPMV)  (46) 

where  pj^  is  the  saturation  pressure  of  methanol  vapor. 

Up  to  this  point,  the  formulation  of  two-phase  mass  transport  in 
all  of  the  regions  of  the  tubular  DMFC  has  been  presented.  To  solve 
the  above  equations,  boundary  conditions  and  additional  constitu¬ 
tive  correlations  are  needed,  which  are  described  in  Tables  1  and  2, 
respectively.  The  cell  geometric  dimensions  and  operating  param¬ 
eters  are  listed  in  Table  3,  and  the  electrochemical  properties  are 
listed  in  Table  4.  A  simulation  code,  which  was  based  on  the  SIM¬ 
PLER  algorithm  with  the  finite-volume  method  [3,25],  was  written 
for  the  above-described  governing  equations  and  boundary  condi¬ 
tions.  It  should  be  noted  that  both  the  main  physical  framework 
and  the  electrochemical  parameters  used  in  the  present  model  are 
based  mainly  on  previous  well-developed  two-phase  models  for 
planar  DMFCs  [3,25].  The  grid  independence  of  the  present  simu¬ 
lation  model  has  also  been  fully  investigated. 

4.  Results  and  discussion 


Fig.  8.  Axial  distributions  of  fluxes  of  (a)  methanol  crossover  and  (b)  water 
crossover  (methanol  solution:  1  M;  current  density:  120  mA cm-2;  anode  flow  rate: 
0.01 5  ms-1;  cathode  flow  rate:  1 .0  m  s_1 ). 


where  the  molar  flux  of  methanol  crossover,  NMLcr,  is  given  by  Eq. 
(23). 

On  the  cathode,  it  is  assumed  that  both  the  cell  current  and 
‘parasitic’  current  are  entirely  consumed  by  the  oxygen  reduction 
reaction  (ORR),  i.e.: 


/  +  /D  = 


!>■ 


jc  dx 


(41) 


where  the  ORR  is  also  given  by  the  Tafel-like  expression: 


jc  —  (1  —  s)/\1<c/gCQ2 


/  Co2  ' 

\  (aJ  \ 

{<%) 

)  p  (  rT^J 

Finally,  the  cell  voltage  can  be  determined  from: 


(42) 


Vcell  =  V0  -  Pa  -  ric  -  / 


(rc  ontact  +  J^) 

\  °mem  J 


(43) 


where  l/0,  ^contact  and  crmem  present  the  thermodynamic  equilib¬ 
rium  voltage  of  the  DMFC,  the  contact  resistance,  and  the  proton 
conductivity  of  the  membrane,  respectively. 

For  the  phase  change  between  liquid  water  and  water  vapor, 
the  rate  of  condensation  and  evaporation  can  be  modeled  using 
the  finite-rate  approach  [11,25]: 

Kvi  =  hvi(ywvPg  -  Pm)  (44) 


4.1.  Detailed  information 

The  distributions  of  several  species,  including  methanol  concen¬ 
tration,  oxygen  concentration,  liquid  saturation,  and  temperature, 
at  a  current  density  of  120  mA  cm-2  for  1  M  solution  are  presented 
in  Figs.  4-7.  The  methanol  solution  and  oxygen  are  fed  at  veloci¬ 
ties  of  0.01 5  and  1.0ms-1,  respectively.  Fig.  4  shows  the  methanol 
concentration  distribution  in  the  anode  porous  region  and  the 
membrane,  as  well  as  the  oxygen  concentration  distribution  in  the 
cathode  porous  region.  It  can  be  seen  that,  due  to  the  electrochem¬ 
ical  consumption  in  the  anode  CL  and  the  methanol  crossover  from 
the  anode  to  the  cathode  through  the  membrane,  the  methanol 
concentration  decreases  sharply  from  the  surface  of  the  anode  DL 
to  the  anode  CL.  While  along  the  flow  direction  in  the  channel,  the 
methanol  concentration  only  shows  a  very  slight  decrease,  which 
indicates  that  the  anode  flow  rate  is  large  enough  to  maintain  a  uni¬ 
form  distribution  of  methanol  along  the  flow  direction.  Due  to  the 
oxygen  consumption  in  the  cathode  CL,  the  oxygen  concentration  is 
decreased  from  the  channel  to  the  cathode  CL,  and  from  the  inlet  to 
the  outlet.  Due  to  the  rather  low  oxygen  transport  resistance  in  the 
gas  phase,  the  decrease  in  oxygen  concentration  from  the  channel 
to  the  CL  is  relatively  small. 

Fig.  5  shows  the  liquid  saturation  distribution  at  both  the  anode 
and  the  cathode  sides.  For  the  anode  side,  the  liquid  saturation  is 
relatively  large  especially  in  the  anode  DL  and  CL  (~0.7-0.95),  while 
in  the  anode  MPL,  the  liquid  saturation  is  rather  small  (~0.2)  due  to 
its  high  hydrophobicity  and  low  permeability.  In  each  layer,  the  liq¬ 
uid  saturation  is  higher  in  the  region  close  to  the  channel  and  lower 
in  the  region  close  to  the  membrane.  That  is  because  carbon  dioxide 
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Fig.  9.  Variations  in  cross-sectional  averaged  (a)  gas  void  fraction  in  the  anode  channel,  (b)  methanol  concentration  in  the  anode  channel,  (c)  oxygen  concentration  in  the 
cathode  channel,  and  (d)  temperature  rise  in  the  cathode  channel  along  the  channel  length  for  various  current  densities. 


gas  generated  in  the  anode  CL  is  transported  from  the  anode  CL  to 
the  channel,  and  is  swept  along  the  channel  by  the  liquid  flow.  Due 
to  the  gas  accumulation,  the  liquid  saturation  is  evidently  decreased 
along  the  channel  direction.  The  accumulation  of  carbon  dioxide  gas 
may  significantly  hinder  the  methanol  supply  from  the  channel  to 
the  anode  CL.  For  the  cathode  side,  the  liquid  saturation  is  rather 
small,  especially  in  the  cathode  MPL  and  DL  (<0.2).  Accordingly, 
the  oxygen  transport  resistance  is  rather  small.  The  liquid  satura¬ 
tion  distribution  and  the  discontinuity  at  the  interfaces  between 
different  porous  layers  in  this  tubular  DMFC  are  similar  to  that  of  a 
planar  DMFC,  as  shown  by  Xu  et  al.  [25].  Also,  contrary  to  the  liquid 
saturation  in  the  anode  side,  it  is  seen  that  there  is  no  evident  dif¬ 
ference  in  the  liquid  saturation  between  the  inlet  and  outlet  region, 
which  means  that  liquid  water  can  be  effectively  swept  out  of  the 
cathode  by  the  air  flow  through  the  cathode  channel.  Previous  find¬ 
ings  can  also  be  clearly  seen  by  plotting  the  radial  distribution  of 
liquid  saturation  at  different  channel  positions  (i.e.,  inlet,  middle, 
and  outlet),  as  shown  in  Fig.  6. 

The  temperature  rise  distribution  across  the  whole  MEA  is  pre¬ 
sented  in  Fig.  7.  As  can  be  seen,  the  temperature  increases  slightly 
along  the  flow  direction  due  to  the  heat  generation  associated  with 
the  electrochemical  reaction.  Moreover,  the  temperature  in  the 
cathode  side  is  generally  larger  than  that  in  the  anode  side.  For 
instance,  the  temperature  rise  in  the  anode  DL  increases  from  0.14 
to  0.9  °C  from  the  inlet  region  to  the  outlet  region,  while  it  increases 
from  about  0.8  to  1.66  °C  in  the  cathode  DL.  That  is  because,  due 


to  the  much  higher  heat  transfer  coefficient  of  liquid  flow  than  gas 
flow,  heat  can  be  removed  more  effectively  through  the  anode  side. 

Fig.  8  presents  the  axial  distributions  of  methanol  and  water 
crossovers  through  the  membrane  from  the  anode  to  the  cathode  at 
the  same  operating  conditions  as  Figs.  4-7.  Both  the  methanol  and 
water  crossovers  are  driven  by  diffusion,  electro-osmosis,  and  pres¬ 
sure  gradient  caused  convection.  As  shown  in  Fig.  8a,  the  methanol 
crossover  is  dominated  by  the  diffusion,  and  the  diffusion  flux 
accounts  for  about  95%  of  the  total  methanol  crossover.  The  convec¬ 
tion  flux  is  negative,  indicating  that  the  pressure  gradient  caused 
convection  is  from  the  cathode  to  the  anode.  Due  to  the  relatively 
low  methanol  concentration  in  the  membrane,  the  convection  flux 
is  rather  small.  Flowever,  it  is  worth  noting  that  the  convection  flux 
may  be  significantly  enhanced  by  adjusting  the  cathode  structure, 
and  thus,  the  total  methanol  crossover  can  be  evidently  lowered.  On 
the  other  hand,  the  methanol  crossover  increases  slightly  and  then 
decreases  slightly  along  the  flow  direction.  That  is  because  along 
the  flow  direction,  the  methanol  concentration  is  decreased,  lead¬ 
ing  to  a  decreased  methanol  crossover,  while  the  temperature  is 
increased,  leading  to  an  increased  diffusivity  and  hence,  methanol 
crossover.  The  trade-off  of  these  two  effects  results  in  the  arch¬ 
shaped  methanol  crossover  distribution. 

Different  from  the  rather  small  flux  of  methanol  crossover,  the 
flux  of  water  crossover,  as  shown  in  Fig.  8b  is  much  larger.  Also, 
although  the  water  crossover  flux  is  dominated  by  both  diffusion 
and  electro-osmosis,  the  effect  of  convection  on  the  total  water 
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Fig.  10.  Variations  in  cross-sectional  averaged  (a)  gas  void  fraction  in  the  anode  channel,  (b)  methanol  concentration  in  the  anode  channel,  and  (c)  temperature  rise  in  the 
anode  channel  along  the  channel  length  for  various  anode  flow  rates  (methanol  solution:  1  M;  current  density:  160  mA cm-2 ;  cathode  flow  rate:  1.0  ms-1). 


crossover  flux  becomes  considerable.  Along  the  flow  direction,  the 
water  crossover  remains  almost  the  same,  which  is  due  to  the  fact 
that  the  species  in  the  CLs,  as  shown  in  Figs.  4-7,  are  almost  the 
same  in  the  axial  direction. 

The  effect  of  current  density  on  the  distributions  of  species, 
including  the  gas  void  fraction,  methanol  concentration,  oxygen 
concentration,  and  temperature  rise,  in  the  channels  are  shown  in 
Fig.  9.  Different  current  densities  ranging  from  30  to  190  mAcm-2 
were  tested.  The  distributions  of  the  cross-sectional  averaged  gas 
void  fraction  in  the  anode  channel  at  different  current  densities  are 
shown  in  Fig.  9a.  It  can  be  seen  that  the  gas  void  fraction  increases 
along  the  flow  direction.  For  instance,  it  increases  from  0%  to  24% 
from  the  inlet  to  the  outlet  for  the  current  density  of  120  mA  cm-2. 
With  the  increase  in  the  current  density,  due  to  the  increased  car¬ 
bon  dioxide  generation,  the  gas  void  fraction  at  the  same  axial 
location  increases  nearly  linearly:  it  increases  from  about  7%  to 
33%  at  the  outlet  when  the  current  density  is  increased  from  30  to 
190  mAcm-2.  Clearly,  the  increased  gas  void  fraction  decreases  the 
area  for  the  mass  transport  of  the  liquid  methanol  solution,  which 
reduces  the  methanol  transport  from  the  channel  to  the  anode  elec¬ 
trode.  The  effect  of  carbon  dioxide  gas  removal  on  the  methanol 
transport  should  be  well  investigated  so  as  to  enhance  the  uni¬ 
form  supply  of  methanol  to  the  anode  CL.  It  should  be  noted  that, 
as  discovered  by  Yang  et  al.  [36],  the  removal  of  carbon  dioxide 
gas  through  the  channel  mainly  forms  bubbly  and  slug  flow.  The 
present  one-dimensional  mixture  model  in  the  anode  channel  can¬ 


not  capture  the  complex  two-phase  flow  behavior,  and  thus,  more 
realistic  two-phase  flow  models  are  needed  to  better  reveal  the 
bubbly  flow  effect  on  the  methanol  transport. 

Fig.  9b  shows  the  methanol  concentration  distributions  in  the 
anode  channel  at  different  current  densities.  Due  to  the  electro¬ 
chemical  consumption  in  the  anode  CL,  the  methanol  concentration 
decreases  along  the  flow  direction.  With  the  increase  in  the  current 
density,  the  methanol  concentration  also  decreases.  Flowever,  the 
decrease  in  the  methanol  concentration  from  the  inlet  to  the  outlet 
is  rather  small,  i.e.,  it  decreases  only  1.3%  at  the  limiting  current 
density  of  190 mAcm-2,  which  indicates  that  the  present  anode 
flow  rate  (0.015  ms-1)  is  large  enough  to  maintain  a  uniform  dis¬ 
tribution  of  methanol  along  the  anode  channel.  Even  at  this  large 
methanol  flow  rate,  the  gas  void  fraction  in  the  channel  outlet  is  still 
considerable,  as  shown  in  Fig.  9a.  Therefore,  the  carbon  dioxide  gas 
bubbles  may  affect  the  methanol  transport  greatly,  and  an  exces¬ 
sive  anode  flow  rate  is  needed  to  avoid  severe  bubble  accumulation 
in  the  channel. 

Similar  to  the  methanol  concentration  distribution  in  the  anode 
channel,  the  oxygen  concentration  in  the  cathode  channel,  as 
shown  in  Fig.  9c,  also  decreases  slightly  along  the  channel  due  to 
the  electrochemical  consumption.  A  higher  current  density  leads  to 
a  larger  decrease  in  oxygen  concentration  from  the  inlet  to  the  out¬ 
let.  Since  the  liquid  saturation  is  almost  zero  in  the  cathode  channel 
(not  shown  here),  and  it  is  relatively  small  in  the  cathode  DL  (as 
shown  in  Figs.  5  and  6),  the  presence  of  liquid  water  in  the  cath- 
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Fig.  11.  Variation  in  the  (a)  cell  voltage,  (b)  flux  of  methanol  and  water  crossover 
through  the  membrane  with  current  density  for  various  anode  flow  rates  (methanol 
solution:  1  M;  cathode  flow  rate:  1.0  ms-1). 


(b) 


ode  has  a  much  lesser  influence  on  the  oxygen  transport  when  the 
oxygen  is  supplied  in  an  active  mode.  The  cross-sectional  averaged 
temperature  in  the  channel  also  depends  on  the  current  density,  as 
can  be  seen  in  Fig.  9d  for  the  cathode  channel.  The  cell  tempera¬ 
ture  increases  along  the  flow  direction  due  to  the  heat  generation 
in  the  MEA,  and  the  increment  becomes  larger  for  higher  current 
densities.  For  the  largest  current  density  of  190  mA cm-2,  the  tem¬ 
perature  in  the  cathode  channel  is  increased  about  2.3  °C  from  the 
inlet  to  the  outlet. 

4.2.  Effect  of  anode  flow  rate 

The  distributions  of  the  cross-sectional  averaged  methanol  con¬ 
centration,  gas  void  fraction,  and  temperature  rise  in  the  anode 
channel  for  various  anode  flow  rates  at  a  current  density  of 
160  mA  cm-2  are  presented  in  Fig.  10.  It  can  be  seen  that  the 
methanol  concentration  (Fig.  10a)  decreases  with  the  decrease  in 
the  anode  flow  rate.  When  the  anode  flow  rate  is  decreased  from 
0.05  to  0.005 ms-1,  the  decrease  in  the  methanol  concentration 
from  the  inlet  to  the  outlet  is  increased  from  about  0.2%  to  3.4%. 
Even  for  the  smallest  anode  flow  rate,  the  change  of  methanol 
concentration  is  still  rather  small.  On  the  contrary,  the  gas  void 
fraction  (Fig.  10b)  clearly  increases  with  the  decrease  in  the  anode 
flow  rate.  When  the  anode  flow  rate  is  decreased  from  0.05  to 
0.005 ms-1,  the  gas  void  fraction  in  the  anode  channel  outlet  is 
significantly  increased  from  about  11%  to  55%.  That  is  because  less 
carbon  dioxide  gas  leaves  the  anode  channel  during  a  lower  flow 
rate.  Therefore,  a  large  anode  flow  rate,  which  may  seem  excessive 


Fig.  12.  Variations  in  cross-sectional  averaged  (a)  oxygen  concentration  in  the  cath¬ 
ode  channel,  (b)  temperature  rise  in  the  cathode  channel  along  the  channel  length  for 
various  cathode  flow  rates  (methanol  solution:  1  M;  current  density:  160mA cm-2; 
anode  flow  rate:  0.01 5  m  s-1 ). 


for  a  uniform  supply  of  methanol  along  the  channel,  is  still  neces¬ 
sary  to  avoid  severe  gas  accumulation  in  the  anode  channel  which 
may  significantly  block  the  methanol  supply  to  the  anode  CL.  For 
the  temperature  rise  in  the  anode  channel  as  shown  in  Fig.  10c, 
due  to  the  heat  removal  by  the  liquid  flow,  the  temperature  is 
increased  along  the  flow  direction.  With  the  decrease  in  the  anode 
flow  rate,  the  heat  removal  per  unit  mass  flow  is  increased,  and 
thus,  the  temperature  is  further  increased.  When  the  anode  flow 
rate  is  decreased  from  0.05  to  0.005  ms-1,  the  temperature  rise  at 
the  channel  outlet  significantly  increases  from  about  0.7  to  7.5  °C. 
Therefore,  a  too  low  anode  flow  rate  not  only  results  in  a  severe  gas 
accumulation  in  the  channel,  but  also  leads  to  a  highly  non-uniform 
temperature  distribution  along  the  flow  direction. 

The  cell  performance,  as  well  as  the  flux  of  methanol  and 
water  crossovers,  for  various  anode  flow  rates  ranging  from  0.05 
to  0.005 ms-1  is  shown  in  Fig.  11.  It  is  found  that  the  anode  flow 
rate  can  greatly  affect  the  cell  performance.  The  cell  limiting  current 
density  increases  with  the  flow  rate,  which  is  similar  to  the  results 
of  a  planar  DMFC  [11].  For  instance,  when  the  anode  flow  rate  is 
decreased  from  0.05  to  0.005 ms-1,  the  limiting  current  density 
decreases  by  about  14%.  Since  the  methanol  concentration  in  the 
channel  (Fig.  10a)  only  decreases  by  about  3%,  the  marked  decrease 
in  the  limiting  current  density  is  mainly  caused  by  the  significant 
increase  in  the  gas  void  fraction  in  the  channel  (Fig.  10b),  which 
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Fig.  13.  Variation  in  the  (a)  cell  voltage,  (b)  flux  of  methanol  and  water  crossover 
through  the  membrane  with  current  density  for  various  cathode  flow  rates 
(methanol  solution:  1  M;  anode  flowrate:  0.01 5 ms-1). 

hinders  the  methanol  transport  from  the  channel  to  the  anode 
DL.  It  should  be  noted  that  investigating  the  effect  of  gas  bubbles 
in  the  channel  on  the  methanol  transport  using  the  present  one¬ 
dimensional  drift-flux  model  is  just  an  initial  approximation,  and 
more  realistic  investigations  need  extremely  sophisticated  two- 
phase  flow  modeling  in  the  porous  region  and  the  channel. 

It  is  seen  in  Fig.  lib  that  the  methanol  crossover  decreases 
sharply  with  the  current  density,  and  approaches  zero  at  the  lim¬ 
iting  current  density.  That  is  because,  with  the  feeding  of  1  M 
solution,  the  methanol  concentration  in  the  anode  CL  significantly 
decreases  with  the  current  density,  which  decreases  the  methanol 
diffusion  through  the  membrane.  At  low  current  densities,  the 
methanol  crossover  is  almost  the  same  for  the  different  anode 
flow  rates,  while  at  large  current  densities,  the  methanol  crossover 
slightly  decreases  with  the  decrease  in  the  anode  flow  rates.  That 
is  also  due  to  the  blockage  effect  of  gas  accumulation  in  the 
anode  channel  on  the  methanol  transport.  Also,  it  is  seen  that  the 
water  crossover  increases  with  the  current  density.  For  instance, 
it  increases  from  about  3.6  to  5.7  p,mol  cm-2  s-1  when  the  current 
density  increases  from  0  to  180  mA cm-2.  This  trend  is  also  similar 
to  that  of  a  planar  DMFC,  which  has  been  well  addressed  elsewhere 
[25].  On  the  other  hand,  since  the  anode  flow  rates  cannot  evidently 
alter  the  water  status  in  the  MEA,  the  water  crossover  is  almost  the 
same  for  the  different  anode  flow  rates. 


cathode  flow  rates  at  a  current  density  of  1 60  mA  cm-2  is  presented 
in  Fig.  12.  From  Fig.  12a,  the  oxygen  concentration  decreases  along 
the  flow  direction  due  to  the  electrochemical  consumption,  and 
the  decrease  becomes  larger  with  the  decrease  in  the  oxygen  flow 
42  rate.  When  the  oxygen  flow  rate  is  decreased  from  1 0  to  0.2  m  s_1 , 

^  the  decrease  in  the  oxygen  concentration  from  the  inlet  to  the  out- 

£  let  is  increased  from  about  0.1%  to  11%.  For  the  temperature  rise 

*|  in  the  cathode  channel  (Fig.  12b),  due  to  the  heat  removal  by  the 

-§  gas  flow,  the  temperature  also  increases  along  the  flow  direction. 

&  With  the  decrease  in  the  oxygen  flow  rate,  the  temperature  is  fur- 

©  ther  increased.  When  the  oxygen  flow  rate  is  decreased  from  1 0  to 

0.2  m  s-1 ,  the  temperature  rise  at  the  channel  outlet  increased  from 
about  0.6  to  2.4  °C.  Therefore,  a  too  low  oxygen  flow  rate  also  leads 
to  a  highly  non-uniform  distribution  of  temperature  along  the  flow 
direction. 

On  the  other  hand,  it  should  be  noted  that  due  to  the  homoge¬ 
neous  model  used  in  the  cathode  channel,  as  well  as  the  Leverett 
relation  used  in  the  cathode  porous  region,  the  predicted  liquid 
saturation  in  the  channel  (<0.01 )  and  in  the  cathode  porous  region 
(<0-2)  is  rather  small.  Accordingly,  the  effect  of  liquid  water  accu- 
9  mulation  in  the  cathode  on  the  oxygen  transport  is  underestimated. 

|  Moreover,  the  change  of  liquid  saturation  with  oxygen  flow  rate  is 
^  also  rather  small,  leading  to  a  negligible  influence  on  the  cell  per- 
«  formance,  as  shown  in  Fig.  1 3a.  It  is  seen  that  the  cell  performance  is 
§  almost  the  same  for  the  different  oxygen  flow  rates  ranging  from  0.2 

2  to  1 0  m  s_1 .  More  accurate  modeling  depends  on  the  measurement 
of  the  capillary  property  of  the  cathode  porous  region,  as  well  as 
the  micro-scale  modeling  of  the  two-phase  flow.  For  the  methanol 
and  water  crossovers  as  shown  in  Fig.  13b,  it  is  seen  that  the  fluxes 
of  methanol  and  water  crossovers  are  all  almost  independent  of  the 
oxygen  flow  rate.  That  is  because,  on  the  one  hand,  the  methanol 
crossover  mainly  depends  on  the  state  of  the  anode  side,  which  is 
almost  unchanged  during  the  increase  in  the  oxygen  flow  rate.  On 
the  other  hand,  the  water  saturation  in  the  cathode  porous  region 
does  not  evidently  change  with  the  oxygen  flow  rate,  which  leads 
to  an  unchanged  water  crossover. 


5.  Conclusions 

Through  a  simple  analysis,  it  is  shown  that  a  tubular  DMFC  sys¬ 
tem  can  yield  a  higher  energy  density  than  a  planar  DMFC  system  by 
properly  choosing  the  tube  radius  and  the  distance  between  adja¬ 
cent  cells.  Furthermore,  based  on  a  steady  state,  two-dimensional, 
two-phase,  non-isothermal  model,  the  mass  and  heat  transport 
characteristics  in  a  tubular  DMFC,  with  different  current  densi¬ 
ties,  methanol  flow  rates  and  air  flow  rates,  were  numerically 
investigated.  The  results  showed  that  along  the  flow  direction,  the 
temperature  slightly  increases,  while  the  temperature  in  the  cath¬ 
ode  side  is  larger  than  that  in  the  anode  side.  Due  to  the  trade-off 
between  the  decrease  in  methanol  concentration  and  the  increase 
in  temperature,  the  methanol  crossover  first  increases  and  then 
slightly  decreases  along  the  flow  direction.  It  is  also  found  that  to 
avoid  severe  bubble  accumulation  in  the  anode  channel,  a  larger, 
excessive  anode  flow  rate  is  needed  to  maintain  a  uniform  distri¬ 
bution  of  methanol  along  the  anode  channel.  Lowering  the  flow 
rates  of  the  methanol  solution  or  air  leads  to  a  higher  tempera¬ 
ture  increase  along  the  flow  channel,  while  resulting  in  negligible 
changes  in  the  methanol  and  water  crossovers  through  the  mem¬ 
brane. 
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